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Abstract. The stellar-dynamical evolution of bound star clusters dur- 
ing the first few Myr is dominated by binary-binary and binary-star in- 
teractions, the rapid sinking of the most massive stars to the centre of 
the clusters and mass loss from evolving stars. The consequences of these 
processes for the binary and stellar population in clusters, and for the star 
clusters as a whole, are studied by following the evolution over 150 Myr 
of a library of compact cluster models containing up to 10^ stars. 



1. Introduction 

The evolution of the primordial binary population, mass segregation, the ejection 
of massive stars, and the effect of mass loss through evolving stars, are dynam- 
ical processes that shape the stellar content and structure of young clusters. A 
particularly interesting problem is the attempt to reconstruct the initial config- 
uration of young clusters, such as the Orion Nebula Cluster and the Pleiades, 
by using the binary population as a memory agent for the dynamical history 
of a population. Since most stars probably form in clusters of some sort, the 
multiple-star properties of the stellar population in the Galactic field is expected 
to be significantly shaped by the dynamical processes in these. 

To understand the appearance of star clusters of any age, their dynami- 
cal and astrophysical evolution needs to be treated accurately and consistently. 
Direct A-body integration with special mathematical methods to treat accu- 
rately and efficiently multiple close stellar encounters as well as perturbed binary 
and higher-order systems, together with modern stellar evolution algorithms, is 
achieved with the NbodyG programme (Aarseth 1999). Using a specially modi- 
fied version of this code, the evolution of clusters with N < 10^ stars is presented 
here. 

2. Models 

Table |l| lists the cluster models. Each was followed for 150 Myr. A is the 
total number of stars, and A^in the number of primordial binaries. Birth binary 
proportions /b = 1,0.6,0 are used with distributions of binding energies as 
in Taurus-Auriga or the Galactic field (Kroupa, Petr, & McCaughrean 1999, 
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hereinafter KPM), respectively. In all cases the central number density has 
the value observed in the Trapezium Cluster, pc = 10^-^^ stars/pc^, and the 
half-mass-diameter crossing time tci = 0.24 Myr. The half-mass radius and 
three-dimensional velocity dispersion are i?o.5 and o"3d, respectively; the initial 
density distributions being Plummer laws. All results shown here are averages 
of A^run calculations per model. In all cases the KTG93 IMF (Kroupa, Tout, & 
Gilmore 1993) is used for 0.08- 1.0 Mq. For all but the lOOflO models, a power- 
law IMF with a = 2.3 (Salpctcr: a = 2.35) for m > IMq, and with a = 0.3 
for 0.01 < m < 0.08 M0, is employed. The lOOflO models have a = 2.7 for 
m, > 1 Mq, and a = 0.5 for brown dwarfs. The mass range is thus 0.01 — 50 
giving the average stellar mass <m> . Random pairing from the IMF is assumed 
to give the primordial mass-ratio distribution. The resulting median relaxation 
time, trei) assumes the cluster either consists of centre-of mass systems (short 
values) or only single stars. Owing to the rapidly changing number of systems 
the true t^ei lies between these two values. 



Table 1. Cluster models 



model 




A^bin 


Ro.5 


<m> 






N 

'run 








[pc] 


[Me] 


[km/s] 


[Myr] 




8fl0 


800 


400 


0.19 


0.37 


1.6 


0.8-1.4 


10 


8f06 


800 


300 


0.19 


0.37 


1.6 


1.0 


10 


8f00 


800 





0.19 


0.37 


1.6 


1.4 


5 


30fl0 


3000 


1500 


0.30 


0.37 


2.5 


2.4-4.4 


3 


30f06 


3000 


1125 


0.30 


0.37 


2.5 


3.0 


3 


30f00 


3000 





0.30 


0.37 


2.5 


4.4 


3 


lOOflO 


10^ 


5000 


0.45 


0.28 


3.2 


14.5 
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3. Stimulated binary evolution 

Binary systems are ubiquitous in all known stellar populations. In particular, 
open clusters such as the Pleiades may contain a binary proportion up to 70 per 
cent (/tot = 0.70) (Kahler 1999). In the Galactic field, /tot ~ 0.6, with a 
Gaussian distribution of logio-P, where P is the binary-star period in days. In 
sparse nearby star- forming regions /b ~ 1, with a primordial period distribution 
that is consistent with rising to the longest periods, and /tot ft 0.5 even in the 
highly concentrated very young Trapezium Cluster (see KPM and references 
therein). Realistic star-cluster models should therefore contain /b > 0.5 at 
birth. 

The evolution of binary orbital parameters through encounters with other 
systems is termed stimulated evolution. During the initial few crossing times the 
weakly bound binary systems are disrupted which can cool the cluster (KPM). 
This occurs for those binaries that have a binding energy, |eb|, smaller than the 
average kinetic energy of the cluster stars, Ck. Binaries with |ebl > ek, on the 
other hand, arc hard, and on average they harden further through interactions 
with cluster stars thereby heating the cluster (Heggie 1975). At any time, bi- 
naries near the hard/soft boundary, with periods P ^ Pth denoting the thermal 
period, are most active in the energy exchange between the cluster field and the 
binary population. The cluster expands as a result of binary heating and mass 
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Figure 1. Illustration of the stimulated evolution of the distribution 
of binary star periods in a cluster (/P = log]^o-f)- ^ binary has orbital 
period Pth when cjsd equals its circular orbital velocity. 

segregation (see below), and the hard/soft boundary, Pthi shifts to longer peri- 
ods (Fig. H). Meanwhile, binaries with P > Pth continue to be disrupted while 
Pth keeps shifting to longer periods. This process ends when Pth^^cut, which 
is the cutoff or maximum period in the surviving period distribution. At this 
critical time, tt, further cluster expansion can be reduced since the population 
of heating sources, i.e. the binaries with P ~ Pth) is significantly reduced. The 
details strongly depend on the initial value of Pth by determining the amount of 
binding energy in soft binaries which can cool the cluster if significant enough. 

Stimulated evolution also leads to a depletion of the mass-ratio distribu- 
tion at small values, but is not efficient enough to thermalize the eccentricity 
distribution, requiring binaries to form with the observed thermal eccentricity 
distribution (Kroupa 1995). 

After the critical time, tt, the expanded cluster has reached a temporary 
state of thermal equilibrium, and further stimulated evolution will occur with a 
significantly reduced rate determined by the velocity dispersion in the cluster, 
the cross section given by the semi-major axis of the binaries, and their number 
density and that of single stars in the cluster. Stimulated evolution during this 
slow phase will usually involve partner exchanges and unstable but also long- 
lived hierarchical systems. The IMF is critically important for this stage, as 
the initial number of massive stars determines the cluster density at t ^ 5 Myr. 
Further binary depletion will occur once the cluster goes into core-collapse and 
the kinetic energy in the core rises. This phase of advanced cluster evolution is, 
however, not addressed here. 

Fig. |2| illustrates the evolution of the binary proportion. All systems are 
counted within and outside the clusters. Initially, the binary proportion decays 
on the crossing time-scale, but slows significantly in all cases after t ~ 1 Myr ~ 
4tcr ~ ^t- The value of f]^ is remembered by 150 Myr for the same IMF. However, 
for model lOOflO /tot crosses over to the /b = 0.6 cases (models 8f06 and 30f06) 



4 


Kroupa 




0.8 


1 1 1 1 1 1 1 1 1 1 

%. 
■ ^CC 


1111 


J 0.6 


8(10 '+ii='^+^^S.|Si:a 




0.4 


, , , , 1 , , , , 1 


30fl0. 

■1 




' ' ' I 

8fl0 



SOflO 



30f06 



1 2 3 50 100 

t [Myr] t [Myr] 

Figure 2. Evolution of the total binary population for models 8fl0 
(thick dashed line), 8f06 (medium dashed line), SOflO (thick solid line), 
30f06 (medium solid line) and lOOflO (dash-dotted line). 
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Figure 3. The number of forbidden orbits at a given time. Models 
are as in Fig. 0. 
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which is due to the smaller number of massive stars. Cluster expansion owing 
to stellar evolutionary mass loss is lessened in model lOOflO, and consequently, 
the clusters retain a higher density. 

Stimulated evolution also leads to the generation of binary systems with 
forbidden orbits. These binaries have an eccentricity, e, that places them above 
or to the left of the upper envelope in e — log^gP diagrams evident in binary 
populations (Kroupa 1995). Forbidden orbits are expected to circularise (i.e. 
eigenevolve) rapidly. This can lead to coalescence in some cases, but when e ~ 1 
after an encounter, the two stars can collide and merge within an orbital period. 
Since the cross section for encounters is much larger for binaries than for single 
stars, this channel leading to stellar collisions will dominate throughout the 
cluster life-time, possibly leading to exotic very massive stars in the early phase 
and/or blue stragglers in older clusters. This is also discussed by Portegies-Zwart 
(2000) for clusters with extraordinarily high central densities and assuming stars 
stick once they touch. 

A more conservative approach is taken here by counting the number of 
forbidden orbits, Oforb, present at discrete time- intervals (Fig. ^). The figure 
shows that larger /b leads to larger Oforb, and that Oforb increases significantly 
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during t^tcr- After this time, the generation of further forbidden orbits is 
reduced, and Oforb is about 0.3 % of the initial number of systems when /b = 0.6, 
but about 0.7 % when = 1. The maximum in Oforb for i~2 Myr when 
N > 3000 is interesting, but remains to be investigated. 



4. Cluster evolution 

Three mechanisms drive internal cluster evolution: (i) energy equipartition (two- 
body relaxation), (ii) binary star activity, and (iii) stellar evolution. Fig. Q 
illustrates the evolution of the central number density, which is defined by the 
number of stars within the central radius of 0.1 pc. The central number density 
decreases on the relaxation time-scale. It also decreases for the clusters with 
/b = 0. This is due to the shedding of binding energy, gained by the massive stars 
as they rapidly sink to the cluster centre as a result of energy equipartition, to the 
cluster field. The central average mass increases by a factor of 2 within 1 Myr, 
and a factor of 4 by 2 Myr, in clusters with the Salpeter IMF above IMq, 
independently of the binary proportion. The cluster (lOOflO) with the steeper 
IMF for massive stars also shows a linear increase of the central average mass, 
but the rate of increase is a factor of 2 slower. 

Returning to the central number density, the clusters with fy, > 0.5 expand 
at a faster rate than the single-star clusters, which is due to heating from the 
significant binary population. The higher density in model 30fl0 (thick solid 
line) than in model 30f06 (medium solid line) for t ^ 0.6 Myr is significant and 
stems from the cooling of the cluster through disrupting wide binaries in 30fl0. 
This effect is seen in reverse in models 8fx (thick and medium dashed lines) 
because the initially smaller a^j) makes more of the binary population effectively 
hard, thus leading to more heating in the case /b = 1. The overall bulk evolution 
is illustrated in Fig. ^. The tidal radius, i?tid) decreases as stars escape, and i?o.5 
increases owing to cluster expansion. The clusters composed of single stars only 
(8f00, 30f00) expand at a slower rate because binary-heating is not effective. 
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Figure 5. Tidal and half-mass radii. Models are as in Fig. ^. 
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Figure 6. The relative number of massive stars (m > 8M0) at dis- 
tances larger than 2 . Models are as in Fig. |^. 

Only the models with N = 800 stars loose approximately 50 % of their mass 
within about 150 Myr. 

Owing to the rapid onset of mass segregation^ stars and systems with lower 
mass have, on average, larger radii than those with higher mass. In particular, 
the calculations show that brown dwarfs have the largest mean radius and are 
preferentially lost from the cluster. However, not only low-mass stars escape. 
The massive stars sink to the cluster core on a time-scale given roughly by 
(<m> /mmass)irei which is ^O.lfrei for the cases studied here. The massive 
stars interact energetically with each other, often leading to their ejection. Many 
calculations show two O and/or B stars leaving a cluster in opposite directions. 
This occurs if at least one of the stars is a binary system that gains binding 
energy and propels itself and the other star out of the cluster. Fig. ^ shows the 
ratio Resc = -^esc/-^now) where A^esc is the number of stars with m > BMq and 
R > 2 iitid , and A^now is the number of such stars remaining in the calculation 
= if iVnow = 0). Thus, between about 5 and 50 Myr, 10-50 % of 
all massive stars are found at distances larger than 2iitid (Fig- This has 
significant implications for studies of the IMF in young clusters, as well as in 
the field. The figure also shows that there is no significant difference between 



ASP Conf. Ser. Style 



7 



the cluster models. However, more massive clusters have a larger Rud, so that 
it takes longer for stars with equal velocity to reach 2 R^^^ for massive clusters. 

5. Conclusions 

Binary-binary and binary-single-star interactions dominate the dynamical evo- 
lution during the first few crossing times, during which time the population of 
soft binaries is depleted significantly, and binaries near the hard/soft boundary 
contribute to cluster heating. The rapid sinking of the most massive stars to the 
centre of the cluster, on a time-scale which is a small fraction of the relaxation 
time, also leads to cluster heating and consequently cluster expansion. About 0.3 
and 0.7 % of the initial number of systems for /b = 0.6 and 1, respectively, are 
hardened sufficiently during an encounter such that efficient circularisation sets 
in. A fraction of such binaries will have their constituent stars collide. Massive 
stars are readily ejected from the cluster cores, and between 10 and 50 % of O 
and B stars are located further than 2 i?tid from their clusters between about 5 
and 50 Myr. The redistribution of binding and kinetic energy leads to cluster 
expansion, and the central density falls on a relaxation time-scale. After about 
4 Myr the massive stars that remain in the cluster start evolving significantly, 
and their combined mass loss leads to further cluster expansion which increases 
the evolution time-scale of the chister. The number of massive stars, and thus 
the IMF, is critical in determining the expansion of the cluster during this phase, 
and by how much the primordial binary population is eroded during the first 
150 Myr. The proportion of primordial binary systems that survive the first 
150 Myr, in clusters that initially have a central number density such as in the 
Orion Nebula Cluster, is 34 and 43 % for = 0.6 and 1, respectively, for a 
Salpeter IMF (a = 2.3) above IMq. If a = 2.7 for m > 1 M©, then /tot ~ 0.35 
despite /b = 1- Finally, the brown dwarf population lies at a larger average 
radius than more massive stars during cluster evolution. 

The results presented here are preliminary in that a larger model library 
is being completed now, and the formidable data reduction is being advanced. 
Many questions remain to be answered. In particular, the effect of gas ejection 
from a newly fledged cluster needs to be addressed. Work on this and other 
issues is in progress. 
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6. Discussion 

David Schaerer: 1) Is there a contradiction between your scenario where mas- 
sive stars are very rapidly ejected from clusters and the scenario of Bonnell, 
Davies, &; Zinnecker (1998) requiring massive stars to form in dense clusters? 
PK: No, not at all. It supports it since the location of massive stars away from 
clusters can be interpreted readily as being due to dynamical ejection. Velocity 
measurements, however, are required to confirm this in any particular case. 
2) Could you comment on the implications of your scenario on the IMF of 
massive stars? Also, Massey finds indications for an IMF less populated by 
massive stars in the field compared to clusters. How does this fit into your 
picture? 

PK: This is a very important question which I am working on. Indeed, the MF 

for massive stars measured in clusters cannot be taken to be the IMF (even if 
stellar evolution is taken account of), and the field stars are likely to have been 
ejected from clusters. Qualitatively this agrees with the finding by Massey that 
the field IMF is steep, because less massive OB stars have a higher chance of 
being ejected with large velocities than the massive ones. However, this requires 
extensive numerical experiments to quantify the ejection velocities and ejected 
mass spectrum, before it can be ascertained that a "primordial field" population 
of OB stars with its own IMF does not exist. 

Hans Zinnecker: I noticed that the time scale for mass segregation you gave 
for the Orion Nebula Cluster was shorter than the crossing time. How can that 
be? If mass segregation is so fast, does that mean that the central Trapezium 
stars may not have been born in the very centre? In other words: do you disagree 
with the recent simulations of Bonnell & Davies (1998)? 

PK: The timescale for mass segregation is not well defined. Spitzer estimated 

the timescale for energy equipartition assuming an isothermal stellar population. 
This estimate is tgqu = {<Tn> /mmass)^rei which can be, formally, shorter than 
tcross- However, tequ is only the timescale for the onset of mass segregation, and as 
soon as the density distribution has changed as a result of energy equipartition, 
the estimate is no longer applicable. Only numerical experiments can lead to 
a better understanding of the mass-segregation time-scale. But little work has 
been done so far on this. Unfortunately Bonnell & Davies used a code that was 
developed by Aarseth for galaxy simulations, and they do not specify the value 
of the softening parameter they adopted. Therefore, it can, at this stage, not be 
said that their calculations correctly treat the cumulative effects of many close 
encounters. One important effect not considered by Bonnell &; Davies is the 
very high primordial binary population. This aids mass segregation because the 
binaries have a significantly larger encounter cross section. 



